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ABSTRACT 

Mestel determined the surface mass distribution of the finite disk for which the circular velocity 
is constant in the disk and found the gravitational field for points in the z = plane. Here we find 
the exact closed form solutions for the potential and the gravitational field of this disk in cylindrical 
coordinates over all the space. The Finite Mestel Disk (FMD) is characterized by a cuspy mass 
distribution in the inner disk region and by an exponential distribution in the outer region of the disk. 
The FMD is quite different from the better known exponential disk or the untruncated Mestel disk 
which, being infinite in extent, are not realistic models of real spiral galaxies. In particular, the FMD 
requires significantly less mass to explain a measured velocity curve. 

Subject headings: celestial mechanics - gravitation - methods: analytical - galaxies: kinematics and 
dynamics 
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1. INTRODUCTION 



iBinnev fc Tremaind ()2008D give a summary of some 
of the explicit potential-density pairs available for the 
study of spiral galaxies and the mathematical methods 
used to discover these relationships. At the best of our 
knowledge, no finite thin-disk solutions are available, and 
this lack has hampered the study of the dynamics of 
spiral galaxies. Solutions on the z = plane exist for 
the exponential disk and fo r the untruncated Mestel disk 
(jBinnev fc Tremainell2008l ) and these have been most fre- 
quently used to model spiral galaxies. A major problem 
is that these disks of infinite extent can over-predict the 
total mass and mass surface density associated with a 
measured veloc ity profile. 

ISchuld (|2009l ) recently presented the complete closed 
form solution of the Maclaurin disk and a few related 
disks over all the space. The Maclaurin disk is applicable 
to the study of spiral galaxies for which the rotational 
velocity rises to the edge of the disk as is seen in irregular 
galaxies. 

Here we present the explicit expressions for the po- 
tential and gravitational field of the finite Mes tel disk 
(herea fter FMD) which was first introduced by IMestell 
jHH). Until now the force vectors of the FMD were 
only known for points in the plane of the disk, and there 
was no analytical solution of the potential. 

Throughout this work we use the notation that M 
is the total disk mass, a is the disk radius, G is the 
gravitational constant, and (i?, z) are the cylindrical co- 
ordinates in each meridional plane. 

2. THE SOLUTION OF THE FINITE MESTEL DISK 

iSchulj ()2009[ ) performed integral transformations of 
the mass distribution and potential of the collapsed 
homeoidal disk with respect to a parameter, determin- 
ing explicit formula for the potential of a family of disks 
which includes the Maclaurin disk. Here we perform a 
similar integration of the collapsed homeoidal disk to find 
the potential and force vectors of the FMD over all the 
space. 

|earIschulz@gmaiI.coin] 



IMestell (|I963f ) constructed the mass distribution of the 
FMD which results in constant circu l ar velo city over the 
entire disk. iLvnden-Bell fc Pineaultl (|1978[ ) gives a con- 
venient form for the surface mass density of the FMD: 



-FMD (it) = -z arecos — 

ZnaK \ a 



R<a. (1) 



Remarkably, the mass distribution given by equation 
([T]) can be obtained by an integral transformation of the 
surface density with respect to a of the homeoidal disk: 
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R<a, 



-'hom 



iR)da. 
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Since this is a linear transformation, the potential and 
gravitational fields of the FMD can be found by applying 
the same transformation to the gravitational field and 
potential of the collapsed homeoidal disk: 
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FMD 



{R,z) 



Fhom{R, z)da, 
^\iom{R, z)da. 



(4) 
(5) 



The explicit expressions for i^hom and ^hom are given in 
ISchulzl Hopl). 

Integration of equation Q can be performed using el- 
ementary functions only. The expressions for the force 
vectors are important because these equations can be 
used to describe the orbits in the gravitational field of 
the FMD. The radial and vertical components of the force 
are: 
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In the z=0 plane equation (|6l) simphfies to: 
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which is ident ical except for sign to the formula given in 
IMestell (Il96l . On the z axis the equation ([7]) simplifies 
to: 

. , GM f\z\\ 

F^,fmd{0, z) = arccot J — ^ . (9) 

az \ a / 

The evaluation of equation ([5]) for the potential of the 
FMD is summarized in the Appen dix. A t variance with 
the force, the potential (equation (|A17p ). cannot be ex- 
pressed using only elementary functions and involves the 
dilogarithm function at complex arg uments. We rec all 
that the Euler dilogarithm function (iMaximonI 120031 ) is 
defined as 

u,iu) = - r '^^dt (10) 



and is a special case of the polylogarithm function defined 
by \jis+i{v) = Jp i^^dt, with Lii{v) = - ln(l - u). 

Figure [1] provides a contour plot of the FMD potential 
at several contour lines which are equally spaced in terms 
of potential difference. The isopotential line which passes 
through the edge of the disk is flattened by a factor of 
about 0.81. 

2.1. Asymptotic Series Expansions 

The potential of the FMD is not analytic on the disk 
and so it is impossible to express the potential as a Taylor 
series expansion. Instead, asymptotic expansions of the 
potential and force vectors were found in terms of power 
series of 1/i?" and 1/z" for i? — > 00 and for 2; — > 00. 

The first few terms of the series expansion of the radial 
force vector in the z = plane are: 

(11) 

Figure [3] shows the behavior of the radial force vector 
near the edge of the disk and also illustrates the conver- 
gence properties of the asymptotic series. As expected, 
the radial force decreases more quickly than Keplerian 
beyond the edge of the disk and quickly approaches the 
asymptotic behavior F oc Including more terms 

increases the accuracy of the partial sum in the region 
near the edge of the disk but convergence is slow and the 
series is not useful for computation. 

The first few terms of the series expansion of the axial 
force vector on the z axis are: 
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Unlike equation pT|) . equation (IT^ is an alternating 
series and has somewhat different convergence behavior. 
The axial force decreases less quickly than Keplerian be- 
yond the limit of convergence at z = a to approach the 
asymptotic behavior F cx 1/z^. 

2.2. Numerical Testing 

Numerical testing was performed in order to verify the 
analytical solution of the FMD. Normalized expressions 
for the potential and force field were defined by choosing 
units such that the disk mass, the gravitational constant, 
and the disk radius were all 1.0. The verification con- 
sisted of the following tests, all of which were performed 
with 80 digit precision: 

• The Laplace equation, 

was verified over a 101 x 101 mesh in the interior 
of the region (0 < i? < 5, < z < 5). Equation 
P3I) was evaluated numerically at each of the mesh 
points using standard finite difference approxima- 
tions and using $ as given by equation (jA17p . The 
derivative of the potential cannot be calculated for 
points i? = or z = and so these points were 
offset by a small amount. The largest calculated 
value of the Laplace equation was 1.9 X 10"^ show- 
ing that equation (fT3|) is satisfied over the entire 
space outside the disk. 

• As shown in the Appendix, the potential of the 
FMD is a real-valued function of real variables 
which involves complex components. The imagi- 
nary parts cancel analytically but, due to roundoff, 
the calculated results might include an imaginary 
component. It was verified that the imaginary com- 
ponents were negligible at all points. 

• The force field is defined by the relationship 

F = — V^. This relation was evaluated numeri- 
cally and the resulting values compared to the force 
vectors given by ^ and ([7]). The evaluation was 
performed at each of the points of the 101 x 101 
mesh described above and the difference was neg- 
ligibly small. 



• From the Laplace equation for a fiat disk 

lim F,{R,z) = -27rGE(i?). 

2— >o+ 



(14) 



This relation was checked numerically at 100 radial 
points and the difference was negligibly small. 

These numerical tests confirm that the analytical ex- 
pressions for the solution of the FMD presented here are 
correct. 

3. DISCUSSION 

A fundamental problem of the study of galaxy dynam- 
ics is to determine the mass model which best fits obser- 
vational data. The FMD is the an appropriate model of 
a spiral galaxy in that both the mass distribution and 
the velocity curve of this disk match those of a typical 
spiral galaxy. 
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• The flat velocity curve of the FMD is typical of 
spiral galaxies. 

Although there is a lot of variation in galaxy rota- 
tion curves, most are flat o r nearl y so in the outer 
parts of the stellar disk. iSofud ()1996D provides 
high-resolution position-velocity curves of nearby 
spiral galaxies based on CO line emission for the 
inner few kpc and on Hi observations for the outer 
region. They show that, in general, the rotation 
curves of spiral galaxies are nearly flat, even in 
the innermost region where the curve rises steeply 
within a few hundred parsecs. 

The flat rotation curve of spiral galaxies might 
be the expected result of evolution. iSaio fc Yoshiil 
()1990l ) demonstrates that a flat rotation curve and 
exponential mass distribution results from the vis- 
cous contraction of a protogalactic disk with on- 
going star for mation. Working i n the context of 
an infinite disk, iBertin fc Lodatd ()1999D , finds that 
viscous flows in a self-gravitating disk and ongoing 
star formation result in a flat rotation curve. 

• The mass distribution of the FMD approximates 
the observed mass distribution of spiral galaxies. 

The mass distribution of spiral galax ies generally 
follow Freeman's law (|Freemanl |197CI[ ) that lumi- 
nosity profiles of disk galaxies usually decrease 
exponentially beyond the inner region: ~ 
exp"^/''". Very often this exponential disk sur- 
rounds a central bulge and a massive central black 
hole. The FMD is a truncated exponential disk 
with a concentrated central region. Figure [2] plots 
the normalized surface mass distribution of the 
FMD given by equation ([T]). Beyond about 0.2 a 
the distribution approaches E oc exp(— i?/ft.o) so 
closely that the two functions are indistinguishable. 



The total mass of the FMD is: 



M 



2aVc 
wG 



(15) 



Where a is the radius of the disk and Vc is the cir- 
cular velocity which is constant over the disk. In 
comparison, an assumed spheroidal mass distribu- 
tion requires a greater mass, by factor of 7r/2 ~ 1.6, 
to explain the same circular velocity profile. 

The scale length of the exponential fit of the sur- 
face mass distribution shown in figure [2] is about 
one third of the radius of the disk: = a/3.1. 
This applies to the outer region of the disk where 
the surface mass distribution of the FMD is indis- 
tinguishable from that of an exponential disk. In 
principle, this ratio could be confirmed by obser- 
vation, noting that the factor 1/3.1 applies to the 
scale length of the mass distribution not the scale 
length of the luminosity function. 

29% of the total mass of the FMD is within the 
inner region R < 0.2a where the distribution ap- 
proaches S cx 1/R. Identifying this region with 
the bulge of a spiral galaxy, it follows that the ro- 
tational velocity can be used to calculate directly 
both the mass of the bulge and the mass of the 



disk. This is significant because these two regions 
typically have different mass-to-light ratios which 
arc difficult to constrain. 

The FMD is a better representation of a spiral galaxy 
than the simple exponential disk with E ~ ^-R/ho ]-,g_ 
cause the exponential disk in infinite in extent and does 
not have a flat rotation curve. Of course, many stud- 
ies have used the so-called "maximum disk" model con- 
sisting of an exponential disk, a bulge component and 
a dark halo component which can reproduce almost any 
rotation curve. However, the procedure of fitting a ro- 
tation curve using these several components is under- 
constrained. The masses of the disk, halo, and bulge 
are free parameters as are the disk scale length and the 
several free parameters which describe the mass distri- 
butions of the bulge and halo. More than one solution 
is possible, depending on the analyst's choice of the free 
parameters. 

The FMD is often conflated with the more commonly 
us ed untruncate d Mestel disk which was also introduced 
in lMesteil ()1963[ ). However the two disks are quite differ- 
ent. The untruncated Mestel disk has a constant circular 
velocity but has infinite mass, is infinite in extent, and 
the mass distribution E ~ <Jq/R does not approximate 
the distribution seen for physical disks. Despite these 
limitations, the untruncated Mestel disk has the appeal 
of simplicity and is often chosen for various studies. 

A fiat disk model requires less mass than a spheroidal 
mass distribution to explain the same rotation curve. Ta- 
ble m compares three disk models which can be used to 
model disk galaxies and shows that a factor of ~ 0.64, 
is required to correct the result of a spherical model to 
that of a FMD for the case of a flat rotation curve. In 
the case of a rising velocity curve, a correction of ^ 0.42 
is required to correct the result a spherical model to that 
of a Maclaurin disk. These large errors which can arise 
from an improper choice of a disk model are often ne- 
glected but should be kept in mind when interpreting 
rotational velocity data. 

The Maclaurin disk is one of the few finite disks which 
have been completely solved ([Schulz 200(1 ) and which, 
along with the FMD, is a reasonable model of a spi- 
ral galaxy. The Maclaurin disk, also known as the 
Kalnajs disk (|Kalnaisl |1972[ ) results from the collapse 
of a rotating spheroid onto the plane. The Maclaurin 
spheroid has been studied since the time of Newt o n (e.g., 
'ChandrasckhaiJ |1987|; iBinnev fc Tremaind [200l iBertinI 
i2000i . for a historical summary and discussions of the use 
of the Maclaurin spheroid and disk to study the stability 
of rotating bodies). 

The Maclaurin disk is defined by 



E = 



3A/ 
27ra2 



v/l-i?7a2 



(16) 



This distribution is flat at the origin and falls off rela- 
tively slowly toward the edge of the disk. The circular ve- 
locity rises linearly to a maximum at the edge of the disk. 
This flat central surface density and linearly rising veloc- 
ity curve is characteristic of irregular galaxies. The total 
mass of the Maclaurin disk can be found from the veloc- 
ity at the edge of the disk with the formula 'iaV^ / {2>itG) . 
In comparison an assumed spheroidal mass distribution 
would require a greater mass, by factor of 37r/4 ~ 2.4, to 
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explain the same circular velocity profile. The large fac- 
tor which must be applied in the case of a rising velocity 
curve has been neglected and has lead to the prediction 
of large amounts of dark matter in irregular galaxies. 

Dark matter halos are not required to explain the 
rot ation curves of the stellar disks of disk galax- 
ies (iPalunas fc Williamsl [2OOOI: [Berth] [2OOOI: [Creze et al.l 
|1998[ )^ The best evidence for dark matter halos comes 
from the velocity curves of the Hi gas which extends be- 
yond the stellar disks. This gas is often seen to have 
a flat rotation curve far beyond the distance at which 
the gravitational attraction of the disk and bulge should 
be Keplerian. There are a few alternate explanations 
of the velocity of extended Hi disks. First, th e MONO 

i Modified Newtonian Dynamics) hypothesis of iMilgroml 
1981 proposes that Newtonian mechanics break down 
below a threshold acceleration qq- Beyond this threshold 
the acceleration due to gravity is proportional to 1/R 
rather than l IR^. Second, the Magnetic Hypothesis 
(iNelsonI [1981 iBattaner fc Floridol I2OOOI : iBattaner et al.l 
l2002f l proposes that the rotation of the Hi gas in the 
outer part of disk galaxies is governed by magnetism 
rather than gravity. This subject is not yet understood 
very well and FMD model presented may help decide 
among the several the several alternative explanations of 
the kinematics of the outer Hi gas. 

3.1. Applications 

The FMD is a suitable model of the galaxy potential 
and force field of typical spiral galaxy and it is expected 
that the solution presented here can be applied to sev- 
eral current astrophysical problems. Some current ar- 
eas of study which would benefit from a more physically 
meaningful model of the gravitational potential of a disk 
include: 



traction of a disk-dominated galaxy decreases with 
increasing height in a way which explains most of 
the observed effect. 



Malasidze_ 

DcionghcJ 



Todual 



de Zeeuw[ ()1988[ ) investigate 



2OOI1 



and 
or- 
bits in the potential of a galaxy disk especially 
as applicable to the study globular clust ers. The 

Jiotential of the disk is assumed to be the iKuzmiiJ 
1956D potential 



= -GM [R^ + {\ z I +zo)' 



and the main parameters of the disk were deter- 
mined based o n a fi t to the measured data of 
iDauphole et al.l (|1996D . A galaxy model based on 
the FMD would provide more relevant information 
for this sort of study, especially near the edge of 
the stellar disk. 

4. CONCLUSIONS 

We have presented the complete closed form solution 
in cylindrical coordinates for the gravitational potential 
and force field of the FMD. The mass distribution and 
the rotation curve of the FMD match those of typical 
spiral galaxies and so this work provide a useful model 
for the study of the properties of disk galaxies. 

The FMD is defined by two parameters, disk mass and 
disk diameter, in contrast to the commonly used maxi- 
mum disc solution which permits many free parameters. 
This simpler description will allow us to standardize the 
characterization of populations of spiral galaxies which 
will aid in the study of, e.g., relations such as the bary- 
onic TuUy-Fisher relation or the Magnitude-Diameter re- 
lation. 



The "galactic fo untain" model o f gas fiow within 
a disk galaxy ([Shapiro fc Field! [19761 : [BregmanI 
[1980a[ [bl) holds that hot gas is ejected from the 
galactic disk by supernova explosions and part 
of this gas condenses into clouds which fall back 
onto the disk at high velocity. This model is 
widely accepted and some progress has been made 
to develop more details of the transport of gas 
and metals and the resul ting chemical e voluti on of 
the disk. See especiallv [Spitoni et al.l ([2008| ) and 
[Spitoni fc Matteuccil ([201]! " and references therein. 

The rotational velocity of extra-planar gas in spiral 
galaxies is obse rved to decrease hnearly with height 
above the disk ljFraternah et al.l [2005[ : [Randl [1997L 
[20051 ) and this be havior has not yet b een under- 
stood completely. [Barnabe et al.[ ([2006[ ) compares 
two models of the extra-planar gas which attempt 
to describe the decrease of velocity with height: a 
ballistic model and a baroclynic model. The bal- 
listic model did not reproduce the data while the 
barocl ynic model was more s uccessful. More re- 
cently, [Marinacci et al.[ ([20111 ) modeled with par- 
tial success the case in which cold 10^ °K) gas 
is ejected from the disk and interacts with hot 
(- 10^ °K) gas in the halo. 

These studies would benefit from the application of 
the FMD model which predicts that the radial at- 



Maple 11 was used to help solve the complicated in- 
tegrals which were encountered and to perform the high 
precision numeric verification. 
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APPENDIX 

THE POTENTIAL OF THE FINITE MESTEL DISK 
Wc start reporting the potential for the coUapsed homeoidal disk ()Schuldl2009D : 

' a + i z 



*hom(i?,2) = TT- 

2a 



arcsm + arcsm 



R J \ R 



(Al) 



where the disk surface density is given in equation ([2]). 

From the hnearity of equation ([3]), it follows that the potential of the finite Mestel disk is found by applying the 
transformation ([3]) to the potential of the collapsed homeoidal disk; 



1 

<i>FMD{R,z) = - $hom(i?,z)da, 

a J a 



(A2) 



-GM /•« arcsin(ii^) , . ... 

= / !^_R_J ^ ' da. (A3 

2a Jo a ' 

'i'FMD is a well-behaved function of the parameter a for arbitrary non-zero R and z. That is, the potential at a 
spatial point (R,z) not on the disk changes smoothly as the disk radius is varied. In fact, $fmd is an analytic function 
of a for complex a over the entire complex plane except for points on the disk. 

We construct the solution to (|A3p using the solution of the simpler integral 

, >N /■'^ arcsin -t- i C) , , , 

M13X)=J \ -dt, (A4) 

where the integration path for equation (|A4[) is along a line segment on the real axis. Here the parameter e was 
introduced to avoid the pole of the integrand at t = and we take C7^0,0<e</3<1, to ensure that the integrand 
is analytic in a region surrounding the path. 
The solution to equation (jA3p is then: 

GM 

$FMD(i?,^) = liiji [Ma/ R,z/R) + M(^/R,-z/R)]- (A5) 
Make the change of variable w — arcsin(t + iC)' 

MPX)= r^^^p^.d^, (A6) 

J wo sm{w) - I C 

where wi = arcsin(a + i^) and wq = arcsin(e 4- The path of integration of equation (jA6p is taken the straight 
line segment connecting points wi and wq. The transformation is one-to-one near the integration path and so this is 
a conformal mapping from one analytic domain to another. 
Expand the integrand of (jA6p using the identities cos(w) = (e*"" + e^'^)/2 and sin(ti;) = (e* — e^*"')/2: 

iwie'^^'^ + 1) 



Define A = VC^TT + C so that 1/A = V?^TT - C: 

iw{e^''" + l 
(e»» + A)(e»--l/A) 



Expand the integrand: 

0i(/3,C) = YTa^/^__ (^-— ^ + — ^ + — ^yy^^jd-. (A9) 
The integral can be solved by using the formulae: 

e* + c 2c c c ^ ' 
dt = (t- l)e* - rfln(l + eVc) - cLia (-eVc) . (All) 
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Equations (|JTO| and (|JTT1) can be proved by differentiation using the definition of Li2, the dilogarithni function, 
given by equation (|T0| . 
Cohect terms: 



MPX) = -iU2{e''"\) - i Li2(-e'"'/A) -tw''/2 + wln{l + e''"/\) + wln(l - Ae''") 



(A12) 



Evaluat e this exp ression, make the substitution A = yl + C^ + Ci ^-i^d simplify the result using the identity 
fGradshtevn fc Rvzhik„2007i 1.622): 

i ln{t — i \/l — t^) ~ arccos(t). 



The solution to equation (jA6p is 
where 



-*Li2{[v/TTTI^TcF-C + */3] (c + v/T+C^)} 

+ arcsin^ (/3 + i C) + [In (2/J) - n2i ] arcsin (/3 + i C) ■ 
The limit at e — )■ is required. After some simplifications we find that: 

lim[</>2(e,C) + <f>2{<^,-C) = iarcsinh^(C) - ^tt^ . 

£->0 

Use equations (jA14p and (jA16|) to evaluate equation (jASP : 



(A13) 
(AM) 

(A15) 



2a 



03(i?, z) + (j)3{R, -z) - i arcsinh^ (^-^^ + -tt^ 



(A16) 



(A17) 



where 



MR,z) 



iLio 



iLi2 



—I arcsin 
2 




(A18) 



The terms of equation (|A17p become large near the z = plane and near the R = axis so that equation (|A17|) is 
unusable those regions. The potential in these regions was found by finding the limits of equation (jA17p . 
The potential in the z = plane is: 



, ttGAI 

$FMD(i?,0) =— In 

Za 

GM 

$FMD(i?,0) 



R 
2^ 

12 



R<a, 



■ arcsin 



i Li2 



-a^ +ia 
R 



The potential on the z axis is: 



$fmd(0,2:) 



iGM 
2a 



In 



i Lio 



I a 



arcsm 



a\ in 



R 



arcsin 



(A19) 



- a'^ +ia 



R> a. 



- Ll2 7 



(A20) 
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Figure 1. Potential of tiic finite Mcstol disk. The isopotential whicfi passes tlirough the edge of the disk is flattened by a factor of 
0.81. The isopotential through (0.2,0.0) is flattened by a factor of 0.53. The potential is flattened near the disk but quickly becomes quite 
rounded. 
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Radial distance 



Finite Mestel Disk Exponential Fit 



Figure 2. Normalized surface mass distribution of the finite Mestel disk. The circular velocity of this disk is constant for > i? < a. The 
dashed line is an exponential which was fit analytically at i? = o/2. The curves are identical in the disk region beyond R > 0.2a 
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FMD Radial Force in the z=0 plane 

0.0 n 




Radial distance 



exact 


— -1/R^2 
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Figure 3. Radial Force Vector of tlie finite Mestcl disk in the plane of the disk. The radial force the disk edge {a = 1) is greater than the 
Keplerian 1/R^ curve by a factor of 1.6. Beyond the disk, the force decreases quickly and joins the Keplerian curve by about R = 1.3a. 
More terms of the series expansion improve the fit somewhat for points beyond the disk edge. 
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Table 1 
Total Mass of Disk Galaxies 
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